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ABSTRACT 

The Wyoming Survey for Ha, or WySH, is a large-area, ground-based, nar- 
rowband imaging survey for Ha-emitting galaxies over the latter half of the age 
of the Universe. The survey spans several square degrees in a set of fields of low 
Galactic cirrus emission. The observing program focuses on multiple Az ~ 0.02 
epochs from z ~ 0.16 to z ~ 0.81 down to a uniform (continuum+line) luminos- 
ity at each epoch of ~ 10 33 W uncorrected for extinction (3cr for a 3" diameter 
aperture). First results are presented here for 98+208 galaxies observed over ap- 
proximately 2 square degrees at redshifts z ~ 0.16 and 0.24, including preliminary 
luminosity functions at these two epochs. These data clearly show an evolution 
with lookback time in the volume-averaged cosmic star formation rate. Integrals 
of Schechter fits to the extinction-corrected Ha luminosity functions indicate star 
formation rates per co-moving volume of 0.009 and 0.014 h 70 M Q yr" 1 Mpc" 3 at 
z ~ 0.16 and 0.24, respectively. The formal uncertainties in the Schechter fits, 
based on this initial subset of the survey, correspond to uncertainties in the cos- 
mic star formation rate density at the >40% level; the tentative uncertainty due 
to cosmic variance is 25%, estimated from separately carrying out the analysis 
on data from the first two fields with substantial datasets. 

Subject headings: galaxies: evolution — galaxies: luminosity function 

1. Introduction 

How do galaxies evolve? Such a broad query covers a wide variety of current research 
interests. One specific approach to answering this question lies in tracking the volume density 
of the cosmic star formation rate as a function of redshift. Early pioneering efforts by Lilly 
et al. (1996), Madau et al. (1996), and Steidel et al. (1996) forged a scenario in which 
the population of star-forming galaxies evolves out to a redshift of approximately unity, and 
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beyond this redshift the star formation rate density plateaus or slowly declines. This general 
picture has since been confirmed through various imaging and spectroscopy surveys carried 
out at ultraviolet, optical, infrared, and submillimeter wavelengths (e.g., Blain et al. 1999; 
Rowan-Robinson 2001; Tresse et al. 2002; Baldry et al. 2002; Wilson et al. 2002; Mann 
et al. 2002; Glazebrook et al. 2003; Hippelein et al. 2003, Heavens et al. 2004, Gabasch 
et al. 2004, Perez- Gonzalez et al. 2005, Le Floc'h et al. 2005, Thompson et al. 2006, 
Babbedge et al. 2006, Sawicki & Thompson 2006, Ly et al. 2007; Shioya et al. 2007; Westra 
& Jones 2007). The rate of this drop-off provides an important constraint to cosmologists- 
understanding the physical processes involving stellar birth, life, and death, and how these 
processes feed back into the interstellar medium is crucial to our overall understanding of 
galaxy formation and evolution (e.g. Somerville, Primack & Faber 2001; Ascasibar et al. 
2002). Galaxy merger rates, the frequency of starburst galaxies compared to active galactic 
nuclei, the production of heavy metals and dust, and the origin of the cosmic background 
levels (Dole et al. 2006) are also closely linked to the cosmic star formation rate. 

While the decrease in the cosmic star formation rate density since z ~ 1 has been 
empirically parametrized as (1 + z) 7 , the strength of this evolution is not well determined. 
Values for the exponent range from 7 pa 1.5 (e.g. Cowie, Songaila & Barger 1999; Wilson et 
al. 2002) to 7 pa 4 (e.g. Lilly et al. 1996; Blain et al. 1999; Tresse et al. 2002, Perez-Gonzalez 
et al. 2005, Le Floc'h et al. 2005, Babbedge et al. 2006), a discrepancy of a factor of six in 
the change in the star formation rate from z = to z = 1. Constraining galaxy evolution 
through measuring changes in the star formation rate density at even higher redshifts is a 
viable and important goal. However, a deep understanding of galaxy evolution, from the era 
of galaxy formation until now, requires an accurate foundation that is based on observations 
of the most recent epochs. 

The current conflict over the strength of the evolutionary trend can be addressed by 
carrying out a uniformly sensitive survey for star formation in hundreds of galaxies at each 
of several epochs. The Wyoming Survey for Ha, or WySH, is a large-area, uniform survey 
for star-forming galaxies at several different epochs that span the latter half of the age of the 
Universd^- The specific aim of the survey is to utilize narrowband imaging of Ha-emitting 
galaxies to conclusively determine the strength of the evolution in the cosmic star formation 
rate density since z ~ 0.8, to provide a firm, low redshift baseline for studies of star-forming 
sources at higher redshift. The purpose of this paper is to introduce the survey, in particular 
the optical imaging program, and to present initial results from the two most recent epochs 
studied. Though only preliminary results are presented in this contribution, the survey will 



1 Results from the z ~ Ha surveys of Gallego et al. (1995), Gronwall et al. (1999), and Hanish et al. 
(2006) will be used for the star formation rate density at the current epoch to supplement our work. 
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ultimately cover volumes at each epoch large enough to minimize cosmic variance due to 
large scale structure and to robustly quantify any dependencies on environment. 

In this paper we outline the optical portion of the survey in § [2J present initial results 
from z ~ 0.16 and z ~ 0.24 in § [3J and summarize in § HI The cosmology assumed is 
H = 70 h 70 km s" 1 Mpc" 1 , Q m = 0.3, and Q x = 0.7. 



2. The Survey 

The optical portion of the survey is being carried out via large blocks of dark time on the 
Wyoming Infrared Observatory 2.3 m telescope (WIR0jE Unlike many previous efforts which 
piece together the star formation history using disparate surveys and techniques/tracers at 
different epochs, the data are being uniformly obtained over the same fields and to the same 
sensitivity to luminosity over the different epochs. The survey is being carried out over 3.5-4 
square degrees from multiple "blank" fields selected to minimize foreground contamination 
(see § I2.2p . By design, these fields are additionally pre-selected to overlap with the target 
areas of deep infrared and ultraviolet surveys. Not only will this survey provide a statistically 
robust measure of the evolving star formation history of the Universe, it will also enable 
important parallel science to be pursued. For example, combining Ha, ultraviolet, and 
infrared data will help us to understand the extent to which the standard individual star 
formation rate indicators are biased (Wilson et al. 2002), and to study the evolution of 
ultraviolet and optical extinction by interstellar dust over time (e.g., Reddy et al. 2006; 
Burgarella et al. 2007; Iglesias-Paramo et al. 2007; Moore et al., in preparation). In 
addition, a large volume, multi-epoch survey will sample well the cluster, group, and field 
environments, leading to strong constraints on the environmental threshold that provides 
the balance between star formation inducement and truncation that occurs in the vicinity 
of clusters (e.g., Dale et al. 2001; Balogh et al. 2002; Finn et al. 2005). 



2.1. Ha as a Star Formation Tracer 

Astronomers rely on a wide range of observational approaches for quantifying star forma- 
tion. The ultraviolet, infrared, and radio restframe continua all provide useful measures of the 



2 The near-infrared portion of the survey (z~0.81) is being completed as a collaborative effort with J. Lee, 
R. Finn, and D. Eisenstein using the Bok 2.3 m and Mayall 4 m telescopes on Kitt Peak. Results from that 
effort will be presented elsewhere. 
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star formation rate, as do line fluxes in the ultraviolet, optical, infrared, and sub/millimeter 
(Kennicutt 1998; Calzetti et al. 2007; Kennicutt et al. 2007). However, for local Universe 
systems Ha photons, indicative of ionized gaseous regions and hence directly related to the 
massive stellar population, provide the conventional standard by which to gauge star forma- 
tion (Kennicutt 1998). Ha is typically the intrinsically strongest optical emission line and 
it suffers less extinction than other traditional indicators that lie at shorter wavelengths. 
Moreover, Ha fluxes are technically much simpler to obtain than star formation measures 
at longer wavelengths. Thus a feasible and direct route to measuring galaxy evolution, one 
that is well-calibrated by local Universe measurements and is accessible from a 2 m class 
ground-based telescope equipped with a wide-field camera, would be to survey the Universe 
in Ha at several different epochs. 



2.2. Field Selection 

In addition to carrying out Ha observations over multiple epochs, we are observing 
multiple fields (Table [Q Figures [Q13]) . Three separate fields help to dampen the effects of 
cosmic variance (Davis & Huchra 1982; Oliver et al. 2000; Huang et al. 2007). Moreover, 
the spacing of the fields in right ascension and their high declinations enable year-round 
observations from Jelm Mountain, located at a terrestrial latitude of +41° and an elevation 
of 9656 feet (2943 m) above sea level, thus maximizing the number of usable nights. These 
fields have been chosen to minimize foreground contamination by zodiacal and Galactic dust, 
bright stars, and nearby galaxies. Our pointings in the ELAIS-N1 and Lockman Hole fields 
overlap with deep surveys that have been carried out in the infrared and ultraviolet. The 
SWIRE Spitzer Legacy project provides maps of these regions at 3.6, 4.5, 5.8, 8.0, 24, 70, 
and 160 /^m (Lonsdale et al. 2003), while the GALEX Deep Imaging Survey offers 30 ks 
integrations in these regions at 1528 and 2271A, integrations that are much deeper than the 
canonical 0.1 ks integration per pointing of the GALEX all-sky survey (Martin et al. 2005). 
Our third field, dubbed WySH 1, conveniently fills the right ascension gap between the 
Lockman Hole and ELAIS-N1, at a northern declination not covered by any other current 
deep field work. 



2.3. Observations and Data Processing 

Optical observations for the survey utilize the WIRO 2.3 m telescope. WIRO is equipped 
with a 2048x2048 CCD camera at prime focus that has C/523 pixels and hence a 17'.9 field 
of view (Pierce & Nations 2002). While most narrowband imaging surveys rely on broad- 
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band imaging for the estimation of the continuum flux levels, we use wavelength-adjacent 
narrowband filters (~ 60A FWHM) to optimize the continuum subtraction. Note that this 
approach leads to a more efficient observing program, since no time is "wasted" on obtaining 
broadband imaging. The filter transmission profiles for the optical portion of this survey 
are provided in Figures H] & [5j The filter widths and central wavelengths in Table [2] are 
computed via 

A = [ T(X)dX (1) 



_ JXT(X)dX 

x = Tnxjdx (2) 

where T reflects the combined transmittance of the filter and the quantum efficiency of the 
CCD, normalized to peak at unity (see Pascual, Gallego, & Zamorano 2007). Note that 
these values are for the typical observing temperature of 0° C. 

The basic unit of observation is a 300 s frame, though several are taken for each filter 
at each location; the individual 300 s frames are slightly dithered to improve image recon- 
struction. IRAF/CCDPR0C is used to run the images through the standard preliminary data 
processing routines. Bias levels and dark counts are removed, the images are trimmed to 
their usable regions (2048x2037 pixels), and the images are flat-fielded using twilight sky 
flats. Fringing affects our images primarily at wavelengths of 8132A and longer, and all im- 
ages suffer from non-uniform illumination that induce a mild halo effect if not removed. The 
amplitude of fringing is at the 10% level or smaller in terms of sky fluctuation, whereas the 
non-uniform illumination results in 1-2% variations in the sky level across the field. CCDPR0C 
is also used to correct for both of these effects, relying on the science frames themselves for 
fringe removal and on median skyflat images for the non-uniform illumination. 

The data presented here were obtained under generally good observing conditions during 
several observing runs spanning 2005, 2006, and 2007. The seeing at Jelm Mountain is 
typically l'/5 FWHM; the observed stellar PSFs from 300 s frames range from l / . / 2-2'/0 for this 
survey. Flux calibration is obtained via 30 s observations of spectrophotometric standards 
(Bohlin, Dickinson, & Calzetti 2001). Based on numerous repeated measurements of the 
same set of standard stars, the uncertainty in the flux calibration is 5%. In addition to 
frequently monitoring the night sky conditions by eye, the photometric quality for each 
frame is based on the relative fluxes in five foreground stars distributed across the field of 
view. Ninety percent of the 300 s frames used in this presentation were obtained under 
photometric conditions (in which foreground stellar fluxes agree to within 2-3% from frame 
to frame). For the remaining 10% of the frames taken in 'nearly photometric' conditions, 
the median flux correction to account for thin cirrus cloud cover is a factor of 1.07. Such a 
correction is applied to each image in conjunction with a correction to convert all fluxes to 
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their equivalent values for zero airmass. 

The multiple 300 s frames for a given field are aligned and stacked to create images 
with longer effective integrations (Table [2]). The 'stack' is an average frame constructed 
from IRAF/IMCOMBINE with avsigclip rejection and a weighting applied to each image 
that is based on the photometric quality described above. Coordinate solutions are written 
to each stack header via IRAF/XYEQ and IRAF/CCMAP, based on the coordinate solutions 
in Digitized Sky Survey imaging and the pixel coordinates for several matched field stars in 
DSS and WySH frames. Inspection of numerous sources spread across the fields of view of 
several stacks shows that WySH coordinates match DSS coordinates to within an arcsecond. 

2.4. Survey Sensitivity 

The sensitivity of the survey is estimated in two ways, with the first based on the 
fluctuation of sky counts and the second based on the observed distribution of detections (the 
identification of sources is described in § |2.5j) . In the first method, assuming a 3" diameter 
aperture and using the observed 3a sky fluctuations of ~ 6.2 ■ 10~ 20 and ~ 3.7 ■ 10 -20 W m -2 
at 7597/766lA and 8132/8199A, respectively, we find that the 3a limiting sensitivity to 
continuum/continuum+line luminosity in our survey is ~ 5 — 7 ■ 10 32 W. 

The noise characteristics within a fixed circular aperture does not necessarily capture 
the full sense of the survey sensitivity. The seeing PSF is of course not constant from night 
to night, nor do we even utilize fixed aperture photometry for our source extraction (§ 12. 5j) . 
Figure [6] shows the relevant figure for the second method of estimating survey sensitivity. 
In this figure we show the observed distribution of signal-to-noise and luminosity for all 
narrowband detections at the targeted redshifts. As can be seen from the distribution, 
the 3cx sensitivity to luminosity is 6 • 10 32 — 3 ■ 10 33 W. In short, the two techniques give 
consistent values for the survey sensitivity to continuum/continuum+line luminosity. The 
sensitivity to line detections is ~ \pl poorer since we derive line strengths via subtraction of 
two narrowband luminosities, as described in the following section. 

2.5. Source Identification and Flux Extraction 

SExtractor (Bertin & Arnouts 1996) is used to identify sources and to extract their 
fluxes. Extensive simulations were performed to optimize the extraction approach for our 
imaging program; key parameter settings include DETECT_THRESH=1.5 and DETECT_MINAREA=3 
0'.'523 pixels. The FLUX_AUT0 output parameter (with the default PH0T_AUT0PARAMS=[2.5,3.5]) 
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is used for the flux extractions, as it effectively recovers simulated fluxes for the variety of 
source morphologies and sizes typical of this survey. Narrowband filter pairs are employed 
at each epoch, and it is unknown in advance whether a particular source will produce line 
emission appearing in 'Filter A' or 'Filter B' of a given filter pair. Thus, SExtractor is used 
in 'dual' mode twice for each stack. First, sources are identified and fluxes are extracted from 
the Filter A stack, then SExtractor is rerun on the Filter B stack using the same positions 
and apertures that were determined from Filter A imaging. The procedure is then repeated 
in the reverse order: identifications and extractions are based on Filter B imaging, and then 
the Filter A image is processed using Filter B positions and apertures. In this way two goals 
are accomplished: fluxes are consistently extracted from Filter A and Filter B imaging, and 
line emitters with weak continuum levels are not overlooked. 

Traditionally, narrowband imaging surveys rely on broadboand imaging to infer and 
subtract off the underlying stellar continuum. Since two wavelength-adjacent but otherwise 
nearly identical narrowband filters are used at each epoch for this survey, the line emission 
is derived from a simple subtraction of the flux from Filter A and the flux from Filter B 
(after accounting for the slightly different filter widths): 

/(Ha + [Nil]) = /a-/b^ if /a>/b^ (3) 

^B ^B 

/(Ha + [Nil]) = / b -/a^ if /b>/a^ (4) 

^A ZXa 

/(Ha) = 0.8/(Ha + [Nil]) (5) 
L(Ha) = 4vr^/(Ha) (6) 
S/N(Ha) = /(Ha + [NH])/v/e(/ A )2 + e (/ B )2 (7) 

where a correction for [N II] that is typical of star-forming galaxies is employed (Kennicutt 
1992; Jansen et al. 2000). The above expressions include the luminosity distance and 
filter fluxes derived as the filter width times the flux density, e.g., /a = Aa/aa- While 
our dual-narrowband imaging is more accurate than a traditional narrowband+broadband 
observational approach, there still may be a selection effect working against low equivalent 
width sources. Hanish et al. (2006) show that 4.5% of their Ha luminosity density comes from 
galaxies with equivalent widths smaller than lOA. Figure [7] displays the Ha equivalent width 
distributions for the two epochs studied in this work, showing that a significant portion of our 
sources have equivalent widths below lOA. The distributions are qualitatively in agreement 
with that presented in the literature (e.g., Gronwall et al. 2004; Haines et al. 2007), though 
there is evidence for larger equivalent widths at redshift z ~ 0.24 than at z ~ 0.16. The 
median equivalent widths are 8.5A and 23A at redshifts 0.16 and 0.24, respectively. 



The redshift of each line emitter is not known a priori and thus detections from galaxies 
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outside of the target epochs could bias to high levels the Ha detection rate and thus the in- 
ferred star formation rate density at each epoch; other prominent optical emission lines such 
as [O III]AA4959/5007, [O II]A3727, and H/3A4861 could be redshifted into the filter band- 
passes (e.g., Fujita et al. 2003). Photometric redshifts are used to cull the contaminators. 
Multiple photometric redshift catalogs are or soon will be available for the fields pursued in 
this work: Rowan-Robinson et al. (2007) from SWIRE, Csabai et al. (2003) from SDSS, and 
Moore et al. (in preparation) from WIRO observations. The SWIRE survey utilizes a suite of 
optical and infrared fluxes along with a variety of spiral, elliptical, and AGN SED templates 
to determine photometric redshifts. The Sloan survey uses a similar technique with ugriz 
data, and our WIRO-based catalog employs UBVRI imaging in addition to our narrowband 
data. In regions where the coverage overlaps (e.g., Lockman Hole), we have compared the 
impact of using photometric redshifts from the various surveys. Survey products such as 
the parameters and integrals of the Ha luminosity functions are fairly insensitive to which 
survey is utilized for photometric redshifts (Moore et al., in preparation). However, for the 
sake of uniformity SWIRE photometric redshifts are used where available, with SDSS and 
our own photometric redshifts employed to supplement the SWIRE coverage. 

The uncertainty on these photometric redshifts is a function of redshift, flux, the num- 
ber and rest wavelengths of the broadband detections, the SED models, etc. The typical 
photometric redshift used in this work has an uncertainty of ~5%. A source is considered to 
be at one of the targeted redshifts if its photometric redshift falls within one of our filter's 
effective bandpasses, i.e., within A/2 of one our filter's central wavelengths (see § 12.31) . The 
sample of Ha emitters is created after applying a 3a cut based on the Ha signal-to-noise. 
The large number of sources detected in our survey leads to a robust statistical representation 
of the star formation rate at each cosmic epoch, but the reliance on photometric redshifts 
implies a certain level of error in tagging individual sources as residing within the targeted 
redshift ranges. The impact of relying on photometric redshifts will be more fully explored 
in a future contribution that presents a larger fraction of the survey. 



A primary tool for understanding the basic characteristics of a galaxy population is the 
luminosity function (Schechter 1976). The Ha luminosity function is constructed for each 
epoch assuming a Schechter profile: 



3. 



Results 
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where a(z) conveys the shape of the function, L#(z) sets the luminosity scale, and <f>*(z) 
represents the overall normalization. The computation of the i th bin of the luminosity 
function and its uncertainty follows 



*(*, bg = Mhilfl , £MZt log Li)] = v^n^L^ 

A log L A log L 

where V(Zj,Lj) is the comoving volume for the j th galaxy in the summation, | logLj — 
log Li |< | A log L, and the bin width AlogL spans 0.4 dex. The Ha luminosity functions 
for each epoch are presented and discussed after a brief description below of corrections that 
are first applied to the data. 



3.1. Accounting for Observational Biases 

3.1.1. Non-Rectangular Filter Transmission Profiles 

If each filter bandpass were to follow a 'tophat' profile, then the comoving volume 
computed in Equation [9] would be simply a product of the filter width A and the number of 
square degrees surveyed (see Table [2]). However, since the filter bandpasses are not tophats, 
it is more difficult to detect faint galaxies over the full bandpass. Hence, the volume in 
Equation [9] for faint galaxies is somewhat smaller than A. The volume correction factors, 
((^filter, where 

V(z,L) cai = ((z) mteT V(z,L), (10) 

can be as large as 10-20%, but the median volume correction factors are Carter— 0.985, 0.981, 
0.990, and 0.988, respectively for filters centered at 7597A, 766lA, 8132A, and 8199A. 

A second effect is related to the non-rectangular nature of the filter bandpass profiles. 
If spectra were available to provide accurate measures of individual galaxy redshifts, then 
each Ha luminosity could be corrected by the inverse of the normalized filter transmission 
at the appropriate wavelength, T(A) _1 . Since such spectra are not available for this survey, 
galaxies that are detected in the low-transmission wings of the filter profiles are measured to 
have artificially faint Ha luminosities. A statistical correction, gaiter , can be applied to all Ha 
luminosities to remedy this effect. The correction is derived by calculating the average inverse 
normalized transmission for each filter, weighted by the survey's distribution of detectable 
Ha luminosities (assuming a Schechter luminosity function). Monte Carlo simulations of the 
effect indicate ^mter ~ 110 ± 0.08. 
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3.1.2. Extinction 



Fluxes are corrected for foreground Milky Way extinction, an effect on the order of 1% 
for these regions of low Galactic cirrus emission. More importantly, the Ha data are also 
corrected for the internal extinction within the survey galaxies. A standard recipe is to 
assume 1 magnitude of extinction at either V band (a factor of ~2.1 at 6563A) or at the Ha 
wavelength (a factor of ~2.5 at 6563A). Another approach is to employ an extinction that 
varies with galaxy luminosity. Hopkins et al. (2001), for example, suggest a transcendental 
equation for the extinction based on the galaxy star formation rate: 



logSFR(Ha) obs = logSFR(Ha) int - a log 



/31ogSFR(Ha) int + 7 
5 



where a = 2.360, (3 = 0.797, 7 = 3.786, and 5 = 2.86 if a Cardelli et al. (1989) reddening 
curve is assumed (see Ly et al. 2007). This approach is adopted for the WySH survey. 
Table [3] lists the extinction correction factor, 9(L) ext , for the relevant bins in the luminosity 
function. The corrected Ha luminosity and star formation rate are thus expressed as 



SFR(H 



L(Ha) cor = 7] mtcI 6(L) cxt L(B.a) ohs 

= 7.9 • 10~ 35 L(Ha) cor [W] 



n) [M Q yr x j 
e(SFR) 



= SFR 



(12) 

(13) 

<-< / 7/1 + e(/fl)7/i + 0.02 2 + 0.15 2 + 0.05 2 + 0.20^4) 



where the Kennicutt (1998) calibration for the star formation rate is utilized. The uncertainty 
in the star formation rate, e(SFR), includes factors of 2%, 15%, 5%, and 20% to account for 
uncertainties in the filter widths, the correction for [N II] (Equation EJ), the flux calibration, 
and the correction factors described in this section, respectively. 



3.1.3. Sample Incompleteness 

An important correction accounts for the increasingly limited ability to detect sources 
with fainter Ha luminosities. An incompleteness correction, k(z, L)^., is derived for each 
bin of the luminosity function, for each epoch observed, 

$(2,logL) cor = K(z,L)r^$(z, log L) ohB (15) 

Extensive Monte Carlo simulations are employed to quantify the impact of incompleteness. 
For each luminosity bin in the luminosity function, IRAF/mkobjects is used to create several 
mock images with the same noise as the actual survey image stacks described in Section 12.31 
These mock images incorporate typical seeing effects and generate sources with distributions 
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similar to those of actual sources in spatial density, equivalent widths, Ha diameters (Dale 
et al. 1999), and of course the luminosity of the particular luminosity bin being simulated. 
Table [3] provides a compilation of incompleteness factors n(z, L) inc for samples that are se- 
lected with a 3<T cut based on the Ha signal-to- noise. The uncertainties in the incompleteness 
are based on the ranges in incompleteness for a variety of simulation parameters (e.g, 500 
vs. 1500 vs 4500 sources in a simulated image; 1 <EW(A) < 100 vs 10 <EW(A) < 50, 
etc). For brighter luminosities the incompleteness reaches an asymptotic value of ~ 0.95 due 
to the effects of sources overlapping one another. Additional pseudo-empirical simulations, 
whereby artificial sources are added to actual image stacks, indicate consistent results with 
the 'pure' simulations. 



3.1.4- Active Galaxies and Star-Forming Galaxies Lacking Ha Emission 

A main goal of this survey is to derive the volume-averaged star formation rate at several 
epochs, Psfr(z). If Psfr(z) is to be interpreted as being solely due to star-forming galaxies, 
then contributions from AGN galaxies will have to be removed. Brinchmann et al. (2004) 
and Salim et al. (2007), for example, respectively use multi- wavelength data for ~ 150,000 
and ~ 50,000 galaxies to derive a ~4% AGN contribution to Psfr( z )- However, these same 
studies estimate between ~1 and 11% of the cosmic star formation rate stems from galaxies 
lacking obvious optical line emission. These two effects have a counteracting impact on 
Psfr{z), and in light of their uncertainties and relatively small impact on the final outcome, 
we do not invoke any corrections to account for them. 



3.2. Preliminary Luminosity Functions at z ~ 0.16 and 0.24 

Figured] shows the initial luminosity functions at z ~ 0.16 and 0.24 for the WySH survey. 
Open circles indicate the data corrected for all issues described above except incompleteness; 
the filled circles also include corrections for incompleteness. The thick solid lines show the 
Schechter fits for all luminosity bins except L(Ha) = 10 33 ' 2 W at z ~ 0.24. The parameters 
for the displayed fits are listed in the first two rows of Table HJ for comparison with the 
literature (e.g., Table 5 of Ly et al. (2007)), we also include in Table @] fit parameters 
assuming 1 and mags of internal extinction at the Ha wavelength. Tentatively, we observe 
a significantly higher luminosity function amplitude <fi* at z ~ 0.24 than at z ~ 0.16. It 
is difficult to compare the characteristic luminosity between the two epochs since it is 
poorly defined at z ~ 0.16. Error bars in Figure [H] reflect the uncertainty in the luminosity 
function amplitude according to Equation [91 summed in quadrature with the uncertainties in 
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the incompleteness corrections. Also included in Figure [8] are Ha-based luminosity functions 
from the literature, largely consistent with our preliminary results. 

3.3. The Cosmic Star Formation Rate Density at z ~ 0.16 and 0.24 

The volume-averaged cosmic star formation rate can be computed by integrating under 
the fitted Schechter function and multiplying by the Kennicutt (1998) star formation rate 
calibration: 



The larger 0* and brighter at z ~ 0.24 lead to an overall larger cosmic star formation 
rate density at that redshift, by a factor of 1.5 for the canonical fits displayed in Figure [HJ 
Table H] provides the integrated cosmic star formation rate densities. For comparison with 
other extinction-corrected, Ha-based cosmic star formation rate densities, we are comfortably 
within a factor of 2 of the values published by Sullivan et al. (2000) at z ~ 0.16 and Pascual 
et al. (2001), Fujita et al. (2003), Hippelein et al. (2003), Ly et al. (2007), Westra & Jones 
(2007), and Shioya et al. (2007) at z ~ 0.24 (see Ly et al. 2007 for values converted to the 
common cosmology also assumed in this work). 

The formal uncertainties presented in Table H] are based solely on the fitting procedure, 
weighted by the uncertainties in the bin data points as described by Equation [91 Also to 
be considered is the impact of "cosmic variance", fluctuations due to the characteristics 
of the particular volume(s) being probed along a survey's line(s)-of-sight (e.g., clusters, 
voids, etc.). An initial estimate of the WySH cosmic variance can be made by performing a 
separate analysis of the two fields studied in this presentation (Lockman Hole and ELAIS- 
Nl). Preliminarily, the cosmic star formation rate densities in these two fields differ by ~25% 
at each epoch, indicating a modest impact due to cosmic variance. 



We report first results from the two most recent epochs of WySH, the Wyoming Sur- 
vey for Ha. The overall scientific goal of the survey is to accurately calibrate changes 
in the Ha luminosity function via <4 square degrees of narrowband imaging from several 
epochs spanning the latter half of the Universe. From a technical standpoint, our aim is 



psFR^roMeyr^Mpc- 3 ) = 7.9 • 10~ 35 £(W) 
where an analytical expression for the luminosity density is 



(16) 




(17) 



4. Summary 
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to achieve relatively small statistical uncertainties by executing a luminosity-uniform survey 
for hundreds of galaxies at each epoch using narrowband filter pairs for improved contin- 
uum subtraction. Preliminary results from approximately 2 square degrees from two fields 
indicate a clear evolution in the volume-averaged cosmic star formation rate, a decrease by 
a factor of ~1.5 from a redshift of z ~ 0.24 to a redshift of ~0.16. Our value at z ~ 0.16 
is comparable to that found from Ha work at z = by Gallego et al. (1995) and Hanish et 
al. (2006), so more substantial evolution appears to have occurred between redshifts of 0.16 
and 0.24 compared to what has occurred after a redshift of 0.16, but additional data are 
needed to more robustly assess this claim. By separately analyzing these initial data from 
the Lockman Hole and ELAIS-N1 fields, we see a ~25% "cosmic variation" in our cosmic 
star formation rate densities at z ~ 0.16 and 0.24. The estimate of the cosmic variance in 
our survey will be more accurately estimated in future work, after more data are added from 
ELAIS-N1 and Lockman, and data are included from our third field ("WySH 1"). 
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and AST0353760) and the Wyoming NASA Space Grant Consortium (NNG05G165H). This 
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in keeping WIRO operational. We acknowledge contributions from WySH affiliates Travis 
Laurance, Mark Reiser, Janice Lee, Rose Finn, and the many students and science teach- 
ers who helped with the observations. We appreciate the early access to the photometric 
redshifts provided by the SWIRE team, and fruitful discussions with Chun Ly. This re- 
search has made use of the NASA/IPAC Infrared Science Archive, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under contract with NASA. 
IRAF, the Image Reduction and Analysis Facility, has been developed by the National Opti- 
cal Astronomy Observatories and the Space Telescope Science Institute. The Digitized Sky 
Surveys were produced at STScI (NAG W-2166). The images of these surveys are based 
on photographic data obtained using the Oschin Schmidt Telescope on Palomar Mountain 
and the UK Schmidt Telescope. Funding for the Sloan Digital Sky Survey and SDSS-II has 
been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the NSF, 
the U.S. Department of Energy, NASA, the Japanese Monbukagakusho, the Max Planck 
Society, and the Higher Education Funding Council for England. 
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This preprint was prepared with the AAS IATj^X macros v5.2. 

Fig. 1. — The ~18' fields targeted at WIRO in the WySH 1 region are shown as squares 
overlaid on an IRAS 100 /im image. 
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Table 1. Target Fields 



Name 


R.A. & Dec. 


/„(100/mi) 


Observing 




(J2000) 


(MJy sr" 1 ) 


Window 


WySH 1 


00 12 40 +32 49 00 


~1.0 


Jul-Nov 


Lockman Hole 


10 47 10 +58 23 00 


-0.9 


Nov-Jun 


ELAIS Nl 


16 11 10 +55 22 00 


~0.4 


Feb-Sep 



Note. — Units of right ascension are hours, minutes, and seconds, 
and units of declination are degrees, arcminutes, and arcseconds. 



Table 2. Optical Filters and Survey Parameters 



Filter 


Filter 


Redshift 


Luminosity 


Current 


Comoving 


Integration 


Sensitivity 


Sensitivity 


A 


A 


for Ho 


Distance 


Area a 


Volume/ Area 


per pixel 


(3a) 


(3a) 


(A) 


(A) 




(/ifo 1 Mpc) 


(□") 


(hf 3 M P c 3 /n°) 


(sec) 


(10- 20 W/m 2 ) 


(10 32 W) 


7597 


60.97 


0.1575 


753 


2.18 


4745 


1200 


6.3 


4.3 


7661 


60.54 


0.1673 


804 


2.18 


5263 


1200 


6.0 


4.7 


8132 


57.68 


0.2392 


1198 


1.92 


9512 


3600 


i.l 


7.1 


8199 


56.55 


0.2493 


1256 


1.92 


10030 


3600 


3.2 


6.1 


8614 


58.56 


0.3126 


1628 




15260 


6000 


2.8 b 


8.7 b 


8687 


64.10 


0.3237 


1695 




17690 


6000 


2.8 b 


9.5 b 


9155 


57.66 


0.3950 


2140 




21920 


9600 


2.2 b 


12 b 


9233 


58.75 


0.4068 


2216 




23390 


9600 


2.2 b 


13 b 



a The projected ultimate coverage for the entire optical survey is 3.5-4 square degrees. 
b Planned. 

°For a 3" diameter aperture, as explain in § 12.41 



Table 3. Ha Extinction and Luminosity Function Incompleteness 



logL(Ha) 




k(z, -t/)inc 


K,(z, -L)i nc 


(W) 




7597/7661A 


8132/8199A 


33.2 


1.18 


0.13+0.02 


0.05+0.01 


33.6 


1.30 


0.71+0.05 


0.27+0.04 


34.0 


1.84 


0.91+0.05 


0.66+0.07 


34.4 


2.38 


0.93+0.04 


0.89+0.05 


34.8 


2.68 


0.94+0.04 


0.94+0.04 


35.2 


3.17 


0.95+0.03 


0.95+0.03 
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Table 4. Luminosity Function Results 



Redshift 


Extinction 


a 


logL* 


log 0» 


PSFR C 




Factor 




(W) 


(ft^Mpc- 3 ) 


/i7oM0yr _:L Mpc- 3 


0.16 


1.18-3.17 a 


-1.62±0.05 


35.5±1.2 


-3.81±0.81 


0.009±0.006 


0.24 


1.18-3.17 a 


-1.55±0.15 


34.7±0.2 


-2.79±0.23 


0.014±0.007 


0.16 


2.51 (1 mag) 


-1.62±0.05 b 


34.7±0.2 


-2.94±0.19 


0.011±0.007 


0.24 


2.51 (1 mag) 


-1.55±0.05 b 


34.8±0.1 


-2.71±0.12 


0.019±0.008 


0.16 


1.00 (0 mag) 


-1.62±0.05 b 


34.1±0.1 


-2.66±0.13 


0.006±0.002 


0.24 


1.00 (0 mag) 


-1.55±0.05 b 


35.0±0.1 


-2.94±0.08 


0.019±0.007 



a See Tabled 
b Fixed; not fitted. 

c The quoted uncertainties are statistical; cosmic variance is tentatively esti- 
mated at 25%. 



Fig. 2. — The ~18' fields targeted at WIRO in the Lockman Hole region are shown as small 
squares overlaid on an IRAS 100 //m image. The 9.3 square degrees mapped at 3.5, 4.5, 5.8, 
and 8.0 jum by the SWIRE program is outlined in bold, while the 9.2 square degrees mapped 
at 24, 70, and 160 /im by SWIRE is outlined with a thinner line. Fields from the GALEX 
Deep Imaging Survey are shown as circles with ~0?6 radii. 



Fig. 3. — The ~18' fields targeted at WIRO in the ELAIS Nl region are shown as small 
squares overlaid on an IRAS 100 fim image. The 11.1 square degrees mapped at 3.5, 4.5, 
5.8, and 8.0 /im by the SWIRE program is outlined in bold, while the 11.0 square degrees 
mapped at 24, 70, and 160 //m by SWIRE is outlined with a thinner line. Fields from the 
GALEX Deep Imaging Survey are shown as circles with ~0?6 radii. 
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Fig. A.— (a) The 7597A/766lA filter pair, (b) The 8132A/8199A filter pair. The filter 
transmission profiles are plotted as filled circles, and the CCD quantum efficiency is shown 
as a dotted line (and enumerated along the righthand axis). The combination of the two 
effects is shown with solid lines. 




Fig. 5.— (a) The 8614A/8687A filter pair, (b) The 9155A/9233A filter pair. The filter 
transmission profiles are plotted as filled circles, and the CCD quantum efficiency is shown 
as a dotted line (and enumerated along the righthand axis). The combination of the two 
effects is shown with solid lines. 
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Fig. 6. — The distribution of signal-to-noise in the individual narrowband images, for sources 
at the targeted redshifts, as a function of narrowband luminosity. A 3a cut corresponds to 
~ 10 33 W. 
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Fig. 7. — The distribution of Ha equivalent widths at z ~ 0.16 (blue dotted line) and 0.24 
(green solid line). 
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Fig. 8. — The preliminary luminosity functions at z ~0.16 and 0.24, based on 98+208 
galaxies observed over ~2 square degrees. The data without incompleteness corrections 
(§ 13 . 1 . 3[) are displayed as open circles, while those corrected for incompleteness are shown 
as filled circles. Error bars reflect the uncertainty in the luminosity function amplitude 
according to Equation [9} summed in quadrature with the uncertainties in the incompleteness 
corrections. The thick solid lines show the Schechter fits for all luminosity bins except 
L(Hcn) = 10 33 ' 2 W at z ~ 0.24; the parameters for these fits are presented in the first two 
rows of Table HI 
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